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It has been suggested recently, based on subtle field-theoretical considerations, that the electro- 
magnetic response of Weyl semimetals and the closely related Weyl insulators can be characterized 
by an axion term BE • B with space and time dependent axion angle 6(r,t). Here we construct a 
minimal lattice model of the Weyl medium and study its electromagnetic response by a combina- 
tion of analytical and numerical techniques. We confirm the existence of the anomalous Hall effect 
expected on the basis of the field theory treatment. We find, contrary to the latter, that chiral mag- 
netic effect (that is, ground- state charge current induced by the applied magnetic field) is absent in 
both the semimetal and the insulator phase. We elucidate the reasons for this discrepancy. 



When a three-dimensional topological insulator (TI) 
[IH3] undergoes a phase transition into an ordinary band 
insulator, its low-energy electronic spectrum at the crit- 
ical point consists of an odd number of 3D massless 
Dirac points. Such 3D Dirac points have been experi- 
mentally observed in TlB^Si-^Se^ crystals [4] and in 
(Bii_ a Jn ;E )2Se2 films [5]. In the presence of the time re- 
versal (T) and inversion (V) symmetries the Dirac points 
are doubly degenerate and occur at high-symmetry po- 
sitions in the Brillouin zone. When T or V is broken, 
however, each Dirac point can split into a pair of 'Weyl 
points' separated from one another in momentum k or 
energy E, as illustrated in Fig. [I] The resulting Weyl 
semimetal constitutes a new phase of topological quan- 
tum matter [6HT4] with a number of fascinating physical 
properties including protected surface states and unusual 
electromagnetic response. 

The low energy theory of an isolated Weyl point is 
given by the Hamiltonian 



hw(k) = bo + vcr • (k — 6), 



(i) 



where v is the characteristic velocity, cr a vector of the 
Pauli matrices, bo and b denote the shift in energy and 
momentum, respectively. Because all three Pauli matri- 
ces are used up in h\y(k), small perturbations can renor- 
malize the parameters, bo, b and v, but cannot open a 
gap. This explains why Weyl semimetal forms a stable 
phase [6 j. Although the phase has yet to be experimen- 
tally observed there are a number of proposed candidate 
systems, including pyrochlore iridates [7J [8] , TI multilay- 
ers [9HT2] , and magnetically doped TIs [T3l [14] . 

The purpose of this Letter is to address the remarkable 
electromagnetic properties of Weyl semimetals. Accord- 
ing to the recent theoretical work [T5UT8] . the universal 
part of their EM response is described by the topological 
#-term, 



S e = 
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E B, 



(2) 



(using H = c = 1 units) with the 'axion' angle given by 




FIG. 1: Low energy spectra in Dirac and Weyl semimetals. 
a) Doubly degenerate massless Dirac cone at the transition 
from a TI to a band insulator. Weyl semimetals with the 
individual cones shifted in b) momenta and c) energy. Panel 
d) illustrates the Weyl insulator which can arise when the 
excitonic instability gaps out the spectrum indicated in c) . In 
all panels two components of the 3D crystal momentum k are 
shown. 



This unusual response is a consequence of the chiral 
anomaly [T9H2T] . well known in the quantum field the- 
ory of Dirac fermions. The physical manifestations of the 
#-term can be best understood from the associated equa- 
tions of motion, which give rise to the following charge 
density and current response, 
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— { bxE-boB). 



(4) 
(5) 



0(r,t) = 2(b-r-b t). 



(3) 



Eq. Q and the first term in Eq. ([5| encode the anoma- 
lous Hall effect that is expected to occur in a Weyl 
semimetal with broken T fTHTQ] . The second term in Eq. 
([5| describes the 'chiral magnetic effect' [22], whereby a 
ground-state dissipationless current proportional to the 
applied magnetic field B is generated in the bulk of a 
Weyl semimetal with broken V . 

The anomalous Hall effect is known to commonly oc- 
cur in solids with broken time-reversal symmetry. In the 
present case of the Weyl semimetal its origin and magni- 
tude can be understood from simple physical arguments 
[TUTU] applied to the bulk system as well as in the limit 
of decoupled 2D layers [18]. Understanding the chiral 
magnetic effect (CME) in a system with non-zero en- 



2 



ergy shift bo presents a far greater challenge. The issue 
becomes particularly intriguing in the case of a Weyl in- 
sulator, illustrated in Fig. [T]i, which will generically arise 
due to the exciton instability in the presence of repul- 
sive interactions and nested Fermi surfaces. According 
to Ref. [15] CME should persist even when the chemical 
potential resides inside the bulk gap. At the same time, 
standard arguments from the band theory of solids dic- 
tate that filled bands cannot contribute to the electrical 
current [23]. We remark that using a different regular- 
ization scheme for the Weyl fermions Ref. [17] found that 
CME occurs in the semimetal but is absent in the insu- 
lator, while Ref. [18] concluded that it only occurs when 
b 2 — 6q > m £>> where denotes the gap magnitude. 
Semiclassical considerations [24] on the other hand pre- 
dict a vanishing electrical current in the Weyl semimetal 
but non-zero 'valley current' proportional to B. 

CME, if present, could have interesting technological 
applications, as it constitutes a dissipationless ground 
state current, controllable by an external field. Disagree- 
ments between the various field-theory predictions, how- 
ever, raise important questions about the existence of 
CME in Weyl semimetals and insulators. The implied 
contradiction with one of the basic results of the band 
theory calls into question whether the results based on 
the low-energy Dirac-Weyl Hamiltonians are applicable 
to the real solid with electrons properly regularized on 
the lattice. In this Letter we undertake to resolve these 
questions by constructing and analyzing a lattice model 
of a Weyl medium. Using simple physical arguments and 
exact numerical diagonalization, we confirm the existence 
of the anomalous Hall effect as implied by Eqs. ( 4|5 ) when 
b 7^ 0. We find, using the same model with bo ^ 0, that 
CME does not occur in either the Weyl semimetal or 
insulator, in agreement with arguments from the band 
theory of solids which we review in some detail. 

Our starting point is the standard model describing a 
3D TI in the Bi2Se3 family, regularized on a simple cubic 
lattice, defined by the the momentum space Hamiltonian 

Ho(k) = 2Xa z (s x sin k y — s y sin k x ) + 2X z a y sin k z 

+ <r x M h , (6) 

with a and s the Pauli matrices in orbital and spin space, 
respectively, and = e — 2t^ a cosk a . For A, X z > 
and 2t < e < 6t the above model describes a strong 
topological insulator with the Z 2 index (1;000). In the 
following, we shall focus on the vicinity of the phase tran- 
sition to the trivial phase that occurs at e = 6t, via the 
gap closing at k = 0. 

It is easy to see that Weyl semimetal emerges when we 
add the following perturbation to Hq, 



Hi(k) = b o- y s z + b • (-(T x s x ,a x s y ,s z 



(7) 
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FIG. 2: The band structure of the Weyl semimetal lattice 
model, displayed along the path k : (7r,0,7r) — >> (0,0,0) — »• 
(0, 0, 7r) — > (7r, 0, 7r). a) Doubly degenerate 3D Dirac point for 
e = 6t and b z = bo = 0. b) Momentum-shifted Weyl point for 
b = 0.9 and bo — 0. c) Energy-shifted Weyl points for b z — 
and bo = 0.7. d) Weyl insulator with b z = and bo = 0.7 
and the exciton gap modeled by taking e = 5.9t. In all panels 
we take X = X z = 1.0, t = 0.5 and the energy is measured in 
units of A. Red circles mark the location of the Dirac/ Weyl 
points. 



follows, V: cr x H(k)a x = H(-k) and T : s y H*(k)s y = 
H(—k). For simplicity and concreteness we focus on 
the case b = b z z, which yields a pair of Weyl points 
at k = ±(b z /2X z )z. The band structure of H = Hq + H\ 
for various cases of interest is displayed in Fig. [2] 

We now address the anomalous Hall effect by directly 
testing Eq. Q. To this end we consider a rectangular 
sample of the Weyl semimetal with a base of (L x L) 
sites in the x-y plane and periodic boundary conditions, 
infinite along the z-direction. The effect of the applied 
magnetic field is included via the standard Peierls substi- 
tution, t — » texp [27ri/^>o // A - dl], where $0 = hc/e is 
the flux quantum, A is the vector potential and the inte- 
gral is taken along the straight line between sites and 
Tj of the lattice. For B = zB{x, y) we retain the transla- 
tional invariance along the z-direction and the Hamilto- 
nian becomes a matrix of size 16L 2 for each value of k z . 
We find the eigenstates (j) n ^ z (x, y) of H by means of ex- 
act numerical diagonalization and use these to calculate 
the charge density 



p(x,y)=e ^2\^n,k z ( x ^y)\ 2 



(8) 
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Nonzero bo breaks V but respects T while b has the 
opposite effect. The two symmetries are generated as 



Figure [3^i displays p for the magnetic field configu- 
ration B(x,y) = $[S(x - L/4) - S(x + L/4)]5(y), i.e. 
two flux tubes separated by L/2 along the x direction. 
In accord with Eq. Q charge accumulates near the flux 
tubes, although p(x, y) is somewhat broadened compared 
to B(x,y). We expect the total accumulated charge per 
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FIG. 3: a) Charge density 5p(x, y) accumulated in the vicin- 
ity of the flux tubes <J> = 0.01$o in the Weyl semimetal. b) 
Total accumulated charge per layer SQ near one of the flux 
tubes, in units of e/27r for indicated values of b z . Dashed 
lines represent the expectation based on Eq |9j). We use 
A = X z = t = 0.5, e — 3.0, L = 14 and L z — 160 independent 
values of k z . Panels c) and d) show the charge accumulations 
as a function of b z in the presence of non-zero Dirac mass and 
bo. Parameters as above except bo = 0.1,0.2,0.3 in c) and 
e = 3.0, 2.9, 2.8, 2.7 for the curves in d) from left to right. 



layer 5Q to be proportional to the total flux, 
e ( b z \ $ 



SQ 



7T \2X Z J $o' 



(9) 



where we have restored the physical units. Fig. ^jp shows 
that this proportionality holds very accurately when the 
flux through an elementary plaquette is small compared 
to <I>o- [When the flux approaches <&q/2 we no longer 
expect Eq. ^ to hold because of the lattice effects.] 
We have also tested the effect of a non-zero Dirac mass, 
vriB — e — 6£, and non-zero bo on the anomalous Hall ef- 
fect. These terms compete with b z and for m 2 D + 6q > b 2 z 
one expects the Hall effect to disappear [17j [18] . This is 
indeed what we observe in Figs. |3]3,d. 

We now address the chiral magnetic effect, predicted 
when bo ^ 0. We consider the same sample geometry 
as above, but now with uniform field B = zB. In order 
to account for possible contribution of the surface states 
we study systems with both periodic and open boundary 
conditions along x. To find the current response we intro- 
duce a uniform vector potential A z along the z-direction 
(in addition to A x and A y required to encode the applied 
magnetic field). The second-quantized Hamiltonian then 
reads 



U(A Z ) = Y,H^(k z - eA z )c\ z 



(10) 



The current operator is given by 
dU(A z 



Jz 



dA z 



^ dH«P(k z ) t 
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This leads to the current expectation value 
dH{k z ) 



dk z 



<t>n,k z ) n F [e n (k z )], (12) 



where rip indicates the Fermi-Dirac distribution and 
£n(k z ) the energy eigenvalues of H(k z ). 



We have evaluated J z from Eq. ( 12 ) for various system 
sizes, boundary conditions, field strengths and parameter 
values corresponding to energy- and momentum-shifted 
Weyl semimetals and insulators. In all cases we found 
J z = to within the numerical accuracy of our compu- 
tations, typically 6-8 orders of magnitude smaller than 
CME expected on the basis of Eq. 

For an insulator, vanishing of J z comes of course 
as no surprise. At T = and using the fact that 
dk z (0n,fc 2 |0n,O = one can rewrite Eq. (12) as 



Jz 
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dk z de n (k z ) 
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dk z 



(13) 



which vanishes owing to the periodicity of e n {k z ) on the 
Brillouin zone. A brief reflection furthermore reveals that 



the integral in Eq. ( 12 ) evaluates to zero even when par- 
tially filled bands are present as long as the system re- 
mains in its ground state or in thermal equilibrium. This 
reflects the well-known fact that one must establish a 
non-equilibrium distribution of electrons to drive a cur- 
rent in a metal, e.g. by applying an electric field. 

There are several notable cases when filled bands do 
contribute currents. A superconductor can be thought of 
as an insulator for Bogoliubov quasiparticles and yet it 
supports a supercurrent. This occurs because Bogoliubov 
quasiparticles, being coherent superpositions of electrons 
and holes, do not carry a definite charge and consequently 
the current cannot be expressed through Eq. (12). In 
quantum Hall insulators non-zero a xy also implies non- 
vanishing current. In the standard Hall bar geometry, 
used in transport measurements, it is well known that the 
physical current is carried by the gapless edge modes, not 
through the gapped bulk. In the Thouless charge pump 
geometry the current indeed flows through the insulat- 
ing bulk but this requires a time-dependent Hamiltonian 
(the magnetic flux through the cylinder is time depen- 
dent). Our considerations leading to Eq. (12) are only 



valid for time-independent Hamiltonians. Finally, there 



where a, j3 represent all the site, orbital and spin indices. 



are known cases [25] when the transition from Eq. ( 12 ) to 
(13) fails because the Hamiltonian is not self-adjoint on 
the space of functions that includes derivatives of (j). This 
can happen when the Hamiltonian is a differential opera- 
tor but in our case H(k z ) is a finite-size hermitian matrix 
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FIG. 4: a) Chiral current J z as a function of energy offset 
bo for various values of the momentum cutoff A. The dashed 
line indicates the field theory prediction Eq. b) The slope 
dJ z /dbo in units of er\ /2tt as a function of cutoff A. Slope 1.0 
is expected on the basis of Eq. 



ior has been conjectured to arise in certain interacting 
systems [28] it is not clear by what mechanism it would 
occur in the ground state of a non-interacting semimetal. 
Our findings indeed confirm the absence of this behav- 
ior in the Weyl semimetal described by a natural lattice 
Hamiltonian. It remains an open question whether this 
fascinating phenomenon can be realized in another quan- 
tum system. 

The authors are indebted to I. Affleck, A. A. Burkov, 
M.P.A. Fisher, A. Grushin and I.F. Herbut for illumi- 
nating discussions and correspondence. The work was 
supported by NSERC and ClfAR. 



with a smooth dependence on k z , which precludes any 
such exotic possibility. In any case, our numerical cal- 
culations addressed directly Eq. (12) so self- adjoint ness 
cannot possibly be an issue. 

Our considerations conclusively establish that anoma- 
lous Hall effect [7HTD] , quantitatively consistent with the 
prediction of the low-energy continuum theory [T5HT8] , 
occurs in realistic Weyl semimetals defined on the lat- 
tice. The chiral magnetic effect [22], implied by the same 
considerations via Eq. (J3|, however runs afoul of the ba- 
sic results of the band theory and is found to be absent. 
Within the low-energy continuum theory the form of the 
axion angle given Eq. (§ can be expected on the basis 
of the Lorenz invariance. In the real solid this symmetry 
is broken at the lattice scale so there is no fundamental 
reason why this form should hold beyond the low-energy 
approximation. To test this hypothesis we have evalu- 
ated current J ; 
i.e. limiting \k- 



from Eq. ([12]) with a momentum cutoff, 

1 



< A in the sum. As shown in Fig. 
when A <C 7T one indeed obtains CME with a magnitude 
consistent with Eq. (J5|. The current however rapidly 
vanishes as the cutoff approaches the extent of the full 
Brillouin zone. We remark that imposing such a cutoff 
has no significant effect on the Hall effect calculation as 
long as A > \b z /2\ z \ because the entire Hall response 
comes from this region of the momentum space [7HTQ]. 
These considerations thus explain the difference between 
the low-energy and lattice descriptions of Weyl semimet- 
als. 

In closing, we note that if correct, the time depen- 
dence of the axion angle implied by Eq. ([3| would en- 
gender some peculiar consequences. One of them follows 
from the Witten effect [26] [27] whereby a unit magnetic 
monopole inserted into the axion medium carries a po- 
larization charge 5Q = — e(n + 6/2tt) with n integer. In 
the Weyl semimetal with a nonzero energy shift bo this 
charge SQ would grow linearly with time according to 
Eq. ([3]). Although such 'quantum time crystal' behav- 
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